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PACS 75 . 45 . + j - Macroscopic quantum phenomena in magnetic systems 

PACS 75 . 50 . Xx - Molecular magnets 

PACS 76 . 30 . -V - Electron paramagnetic resonance and relaxation 

Abstract. - Spin relaxation between the two lowest-lying spin-states has been studied in the 
S = 4 single molecule magnet Ni4 under steady state conditions of low amplitude and contin- 
uous microwave irradiation. The relaxation rate was determined as a function of temperature 
at two frequencies, 10 and 27.8 GHz, by simultaneously measuring the magnetization and the 
absorbed microwave power. A strong temperature dependence is observed below 1.5 K, which 
is not consistent with a direct singie-spin-phonon relaxation process. The data instead suggest 
that the spin relaxation is dominated by a phonon bottleneck at low temperatures and occurs by 
an Orbach mechanism involving excited spin-levels at higher temperatures. Experimental results 
are compared with detailed calculations of the relaxation rate using the universal density matrix 
equation. 
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Quantum tunneling of magnetization (QTM) [1,2] and 
quantum phase interference [3, 4] have been intensively 
studied in single molecule magnets (SMMs). These ma- 
terials have also been suggested as candidates for qubits 
in quantum processors [5] and for applications in molec- 
ular spintronics [6]. Quantum manipulation of the spin 
in such materials can be envisioned provided they have 
Sufficiently long spin relaxation (Ti) and dephasing times 
{T2). A first step has been the study of coherent QTM, 
in which the tunneling rates are faster than the rate of 
decoherence [7, 8] . These studies provided a lower bound 
on T2 of about 0.5 ns. In SMM crystals dipolar interac- 
tions may limit the dephasing time, making it necessary to 
work with dilute SMM ensembles [9-11] or even individ- 
ual molecules. For instance, a recent study of dilute doped 
antiferromagnetic wheels demonstrated a phase relaxation 
time of several microseconds at 1.8 K [12]. 

The phase relaxation rate r2 — I/T2 is generally much 
greater than the energy relaxation rate F. For instance, 
an upper bound on T2 of 50 ns was found in dilute Ni4 
solutions at 130 GHz and 5.5 K, where energy relaxation, 
rather than dipolar or hyperfine interactions, may be the 
limiting mechanism [10]. Interestingly, the relaxation rate 



F due to a direct spin-phonon process of a SMM embed- 
ded in an elastic medium can be derived without any un- 
known coupling constant [13,14]. Moreover, collective re- 
laxation effects are expected in SMM single crystals, such 
as phonon supcrradiance [15]. 

Direct measurements of the magnetization under pulse 
microwave (MW) irradiation enable study of spin- 
dynamics in SMM crystals [16-19]. In Feg, long pulses 
(> 10 /is) of resonant MW radiation drive spins and 
phonons out of equilibrium, and heating effects or phonon 
bottleneck (PB) preclude the observation of fast dynam- 
ics [16,18]. However, using very short high power MW 
pulses, refs. [17] and [19] showed that these effects can be 
circumvented, enabling the observation of spin-dynamics 
on microsecond time scales. Ref. [17] explains the observed 
temperature dependence of the level lifetimes with direct 
spin-phonon coupling whereas ref. [19] concludes that a 
PB develops and plays an essential role in the magnetiza- 
tion dynamics. 

Several experiments on SMM crystals in the quasi-static 
regime [20,21] were explained based on the ad hoc model 
of PB [22,23]. Only recently has this effect been investi- 
gated on microscopic basis [24,25]. When Ns, the num- 
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ber of spins in a crystal, becomes comparable to A^ph, the 
number of resonant phonon modes (which can take energy 
out from the spin system), spin relaxation leads to the 
so-called bottleneck plateau [24]. Only because phonons 
themselves can relax to another reservoir is thermal equi- 
librium eventually attained. Hence, the collective dynam- 
ics of spins plus resonant phonons depends on phonon re- 
laxation, which can be limited at low temperatures [26]. 
In SMM crystals the bottleneck parameter B = Ns/Nph 
is very large (^ 10), even though the resonant lines are 
broadened by inhomogeneities [25]. Hence, collective ef- 
fects in spin-phonon relaxation should be strong in SMMs. 

This Letter presents a study of magnetization relaxation 
in the SMM Ni4 under steady state conditions of low am- 
plitude and continuous microwave irradiation. The gen- 
eral methods used can be applied to any spin system, in- 
cluding other molecular magnets. A large transverse field 
is used to tune the energy splitting between the two low- 
est spin-states [8] , while small deviations from equilibrium 
magnetization induced by a resonant MW field are mon- 
itored. Simultaneous measurements of the absorbed MW 
power and of the magnetization enable measurement of en- 
ergy relaxation rate as a function of temperature. Exper- 
imental results are compared with detailed calculations of 
the relaxation rate between the first excited state and the 
ground state, using the density matrix equation (DME) 
with the relaxation terms in the universal form [13, 14]. 

In steady state, the energy level populations are con- 
stant. Under continuous wave (cw) irradiation of fre- 
quency / = ti;o/(27r), the absorbed power Pabs is pro- 
portional to the spin-photon transition rate A and to the 
population difference Nq = Na — Ni, between the states 
I a) and \b) involved in the transition. The latter is re- 
duced from its equilibrium value iVg"' by the resonant 
perturbation. This reduction depends on the ratio be- 
tween the transition and relaxation rates. For a two level 
system, the total relaxation rate is F = F^f, -I- Ti,a and 
Pabs = hcooANo = hujQANl^'^/{l + 2A/F). A can be elim- 
inated giving F ~ 2Pabs/(/i'^o(-^o'' — -^o))- The popula- 
tions can be obtained from magnetization measurement, 
M = Nauia + NbTUb = M°i + AM. Thus the total relax- 
ation rate between \a) and |6) can be expressed in terms 
of the measurable quantities Pabs and AM as follows: 
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Eq. (1) is a simple result of the balance between absorp- 
tion and relaxation in steady state for a two level system. 
If the system has more than two levels - with levels a 
and h lowest in energy - one can still consider an effective 
two- level model in which relaxation between \a) and |6) is 
assisted by upper levels, provided that the temperature is 
sufficiently low so that the populations of all other levels 
are small. The use of eq. (1) is still justified because the 
upper levels do not contribute to Pabs and AM. 

Simultaneous measurements of Pabs and AM have been 
previously used to measure relaxation times in paramag- 




Figure 1: (Color online) Schematic of the integrated sensor 
used in this study. 



nets and ferromagnets [27,28]. Recently, quantitative de- 
termination of the magnetization suppression under high 
power cw irradiation was obtained in the SMM Feg us- 
ing a commercial SQUID [29]. Here, we use a novel in- 
tegrated sensor that incorporates a microstrip resonator 
with a micro-Hall magnetometer [30], which enables low 
power studies. Briefly, the Hall sensor is fabricated from a 
2 dimensional electron gas in a GaAs/AlGaAs heterostruc- 
ture. A microstrip resonator is capacitively coupled to two 
transmission lines and lithographically placed on the Hall 
sensor, as shown in fig. 1. A thin dielectric separates the 
resonator from the heterostructure material. The mag- 
netometer, biased with a dc current of 10 /xA, permits 
detection of changes in the magnetization AM as small 
as lO""* of the saturation magnetization Alg of our SMM 
crystals. The resonator has a fundamental resonance fre- 
quency of 10 GHz, for which the ac field /lac is maximum 
at the sample location. This is also the case for the third 
harmonic, at 27.8 GHz. The large filling factor of the res- 
onator enables the detection of absorbed power in the nW 
range for « 100^/im'^ samples, ft-ac was calibrated at room 
temperature at each frequency through DPPH saturation. 
A vector network analyzer is used as a MW source and 
for transmission measurements. The integrated sensor is 
mounted in a '^He cryostat with a base temperature of 
0.37 K, in a superconducting vector magnet that enables 
magnetic fields to be applied in arbitrary directions with 
respect to the axes of the crystal. 

[Ni(hmp)(dmb)Cl]4 (Ni4) is a particularly clean SMM 
with no solvate molecules present in its crystal phase and 
only 1% (natural abundance) of nuclear spins on the tran- 
sition metal sites [31]. This results in narrower EPR peaks 
than in many SMMs [32] . The spin Hamiltonian of Ni4 is 
to first approximation: 

Hs = -DSl - BSt + CiSl + St)- ^ibH gS, (2) 

where the first term is the uniaxial anisotropy, the second 
and third terms are 4*''-order anisotropy terms, and the 
last term is the Zeeman energy. The 5* = 4 ground state 
of the molecule at low temperature is a consequence of 
ferromagnetic exchange interactions between the four Ni^^ 
{S = 1) ions. The uniaxial anisotropy leads to a large 
energy barrier DS^ « 12 K to magnetization reversal be- 
tween states of projection Sz = ±4 along the easy axis 
of the molecule. The large C-term is at the origin of the 
fast tunneling (« 10 MHz) observed at zero field [33] and 
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Figure 2; (Color online) (a) Energy levels of Ni4 vs. longitudi- 
nal field H^ in the presence of a transverse held H± = 3.1 T. 
(b) Absorbed power and (c) magnetization changes vs. H^ 
measured with an the integrated sensor as the Ni4 crystal is 
irradiated at 27.8 GHz. The dashed vertical lines show the 
resonant transitions between the two lowest levels \a) and jb) 
corresponding to this frequency. 



a recent analysis shows that the 4*''-order terms in eq. (2) 
result from a finite ratio of the (second order) single ion 
anisotropics to exchange constant [34] . 

A pyramidal Ni4 crystal containing Ns ~ (7±0.5) x 10^^ 
molecules is affixed to the integrated sensor with vacuum 
grease and oriented so that /lac is within 20° of z, the 
axis of the pyramid. A constant transverse field H± is 
applied^ perpendicularly to z, while the longitudinal field 
Hz is swept [8,30]. In this configuration, an avoided cross- 
ing between the two lowest energy levels \a) and \b) oc- 
curs, with a tunnel splitting A (fig. 2(a)). The energy 
levels and eigenstates are calculated by direct diagonal- 
ization of the Hamiltonian (cq. (2)). The character of 
the states \a) and \b) changes with H^- At H^ = 0, 
they have zero magnetization along z. To have a mea- 
surable magnetization we work at finite Hz, where rria 
and TO(, are also finite and almost opposite. In the pres- 
ence of cw MW irradiation, power absorption can be de- 
tected when the splitting between \a) and |6) matches the 
photon energy fkuo (fig. 2(b)). Simultaneous magnetiza- 
tion changes AM are observed (fig. 2(c))^. By repeat- 
ing this experiment at different transverse fields with the 
MW frequency fixed, it is possible to determine the zero 



^At an angle <p ~ 11° with the hard axis, estimated from the 
crystal geometry [33]. 

^Ni4 molecules in a different molecular environment are respon- 
sible for the peak shoulders at ai ±0.15 T [32]. 
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Figure 3: (Color online) (a) Equilibrium magnetization, (b) ab- 
sorbed power and (c) magnetization changes vs. temperature, 
measured for two values of the applied field, corresponding to 
transition frequencies of 10 and 27.8 GHz. The continuous lines 
are calculations (see text). 



field splitting parameters [8,30], which are found to be in 
good agreement with those deduced from high frequency 
EPR experiments [32,33]. Using the two frequencies avail- 
able for our microstrip resonator, we obtain D = 0.75 K, 
B = 7 X 10^3 K, C = 2.9 X lO""* K, g^ = g.y = 2.23 and 
gz = 2.3. 

We now concentrate on the behavior of Pabs and AM 
vs. temperature for two values of the applied field (fig. 3). 
{H±, Hz)i = (2.55, 0.03) T corresponds to a transition fre- 
quency of 10 GHz and {Hj_,Hz)2 == (3.1,0.11) T to 27.8 
GHz. In both cases, only the two lowest-lying states are 
significantly populated within the temperature range in- 
vestigated. At 1.5 K, less than 3% of the molecules occupy 
other states. The first such state lies several Kelvins above 
b (see fig. 2(a)), and thus the effective two-state model dis- 
cussed above is valid. The temperature dependence of the 
equilibrium magnetization in the absence of MW, M'^'^, is 
plotted in fig. 3(a). It is well reproduced by the two-state 
model (continuous lines): 



~m7^ S{l + cxp{-hiJo/{kBT))) ' 



(3) 



where rua, mt, and ujo are calculated from diagonalization 
of Hs. It is found that ma/S = 0.45 and mb/S = -0.432 
for {H±,Hz)i, and that rUa/S = 0.425 and mt/S = 
-0.329 ior{H^,Hz)2. 

When irradiating the sample with MW, the power was 
deliberatly kept small in order to investigate small devia- 
tions from equilibrium, typically less than 1%. The lowest 
achievable sample temperature under these conditions is 
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0.385 K, only 15 mK above the cryostat base temperature. 
The measured temperature dependence of Pabs is plotted 
in fig. 3(b) together with the formula: 
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Fj is an upper bound on the dephasing rate, and can be ex- 
tracted from the line width, found to be about 4 GHz [30] . 
The only fitting parameter for the continuous lines dis- 
played in fig. 3(b) is /lac, which is found to be 1.9 mOe at 
10 GHz and 1.3 mOe at 27.8 GHz. These values are 30 to 
50% less than expected from calibration. The difference in 
the temperature dependences of Af'^ and Pabs under the 
different applied fields is simply due to the different level 
splittings, ftwo, in the Boltzmann factors of eqs. (3) and 
(4). Finally AM. the magnetization change induced by 
the MW absorption, is plotted in fig. 3(c) for the two ap- 
plied fields. Its temperature dependence is much stronger 
than that of M'^'^ and Pabs- ^M/Mg decreases by almost 
two orders of magnitude at 10 GHz as T increases from 
0.385 to 1.2 K, and by almost a factor 20 at 27.8 GHz. 

With the data presented in fig. 3 and eq. (1) we have 
extracted^ the energy relaxation rate F. It is displayed in 
a semilog plot vs. the inverse temperature \/T in fig. 4. 
The general trend of the data at 10 and 27.8 GHz is a 
moderate increase of F with temperature at low temper- 
atures (T < 0.8 K) followed by an exponential growth at 
higher temperatures. 

A direct spin-phonon relaxation process from \b) to \a) 
is not consistent with the data, because it would lead to 
a more gradual temperature dependence of the relaxation 
rate. In the case where the transverse field is mainly re- 
sponsible for the energy splitting between states, it reads 
[14], 
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where Q ^ 1 when SH, < Hj_ and Et = {pvlh^y/'^ is 
the characteristic energy of the spin-phonon interaction. 
p is the density of Ni4, and Vt the speed of the trans- 
verse phonons. The temperature dependence of Fsp,T is 
much weaker than that of the measured rate F, even at 
low temperature. The zero temperature spin-phonon re- 
laxation rate Fgp given by the prefactor of cq. (5) is very 
sensitive to the precise value of the transverse speed of 
sound, since vt appears to the fifth power. The latter was 
extracted from the heat capacity C{T) of Ni4 [35] using 
the extended Debyc model [36]. The soft phonon mode, 
Vt = 940 m/s, is nearly transverse and the most impor- 
tant for the spin-phonon relaxation, whereas the two stiff 
modes with v = 2250 m/s can be neglected. The prefac- 
tor i in eq. (5) accounts for the fact that only one phonon 
mode contributes to the spin-phonon process. Using the 



^In order to extract T up to T = 1.6 K, the 10 GHz data of fig. 4 
also include measurements at a MW power of Pin = 1 /^W. 
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Figure 4: (Color online) Dependence of the energy relaxation 
rate at 10 and 27.8 GHz on inverse temperature deduced from 
the simultaneous measurements of fig. 3 and from eq. (1). The 
continuous lines are calculations using the universal DME and 
the dashed lines are fits of the data including the Orbach mech- 
anism and PB (see text). 



density of Ni4, p = 1.4 g/cm^ [37], we obtain Et = 75.7 K. 
Using this value in eq. (5), the expected zero temperature 
relaxation rate from the direct process is found to be 40 
and 200 times larger than the observed rate at 10 GHz and 
27.8 GHz, respectively. Such a large disagreement can not 
be assigned to an uncertainty of Vt, which would have to 
differ by a factor 2 to 3 from the value we have taken. Also 
processes other than rotation of magnetic molecules with- 
out distortion, unaccounted for in the universal form of 
the spin-phonon coupling, cannot explain the discrepancy 
since they can only increase the relaxation rate. 

This suggests that a PB is necessary to explain the data. 
In fact, an excellent agreement with the data at each fre- 
quency is obtained with the form of relaxation rate that 
combines PB and Orbach mechanism (see dashed lines in 
fig. 4): 
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This temperature dependence of the relaxation rate has 
also been observed in dilute paramagnetic salts [38] . 

The temperature dependence of the PB mechanism in 
the presence of phonon relaxation Fph can be obtained 
from eqs. (60) of ref. [24] by including a large inhomoge- 
neous broadening, as done in ref. [25]. In that case, the 
bottleneck parameter has the form: 



N pph(t^o) 



Ns 2TT''n% 



(7) 



Ny/2TT i^qSujo ' 
Here, the spin density of states ps{uJo) is assumed to have a 
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Gaussian line shape of width Sloq, N is the number of unit- 
cells and vq their volume, ilu — Vi/vq is proportional 
to the Debye frequency. Integrating out rt and then over 
the narrow natural line shape in eqs. (60) of ref. [24] with 
the added inhomogeneous broadening leads to equations: 
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that are similar to those of refs [22,23]. p^g is the pop- 
ulation of the spin's excited state and n^^g the phonon 
population. From linearization of eqs. (8) for Bi^g ^ 1, it 
can be found that the lowest rate governing the relaxation 

is: 

rpiJsp.T 
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Bujo,T^sp,T + Tph 
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where B^^^t = -Bi^j, tanh {hujo/{2kBT)). In the strong 
bottleneck limit, Fph ^ B^g^T^sp,T, the so-called bottle- 
neck rate A_ ^ Fb = Fph/-Bwo,T dominates, as in our 
experiments at low temperature, i.e., the first term of 
eq. (6). In Ni4, the unit-cell contains four molecules and 
its volume is vq = 5777.1 A^ [37]. Hence, Ns/N = 4 



and no = 8.3 x 10" s-\ 

B^„ w 180000 at 10 GHz and B^,, 



Using SiOo = 25 X 10^ s 
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24000 at 27.8 GHz, 
values which also depend strongly on vt . From the prefac- 
tors of the PB mechanism used to fit the data (see dashed 
lines on fig. 4), we find that in our experiments, the phonon 
relaxation time Tph — 1/Fph ~ 0.4 iis at 10 GHz and 
« 0.8 fis at 27.8 GHz. This corresponds to a phonon mean 
free path of about 500 /im, for a crystal whose largest 
dimension is 180 /xm. These values are comparable to 
published ones [23], and reveal that the phonons undergo 
several reflections on the crystal walls before escaping. 

To better understand the relaxation including the Or- 
bach mechanism that contributes to cq. (6), we have per- 
formed calculations of the relaxation rate between the two 
lowest-lying states in the whole temperature range, using 
the density matrix equation (DME) with the relaxation 
terms in the universal form [13, 14]. In this approach, the 
complete set of energy levels and the effect of temperature 
are taken into account. The explicit form of the crystal 
fields enters through the eigenfunctions and eigenvalues 
of the Hamiltonian (eq. (2)) that are universal building 
blocks of the relaxation term of the DME. This makes in- 
clusion of various anisotropics easy since it does not change 
the structure of the latter. The continuous lines in fig. 4 
are the results of DME calculations for the applied fields 
{H±,Hz)i and {H±,Hz)2 used in the experiments. The 
Hamiltonian parameters and characteristic energy of the 
spin-phonon coupling needed for the calculation arc the 
same as given above. The low temperature part of the 
DME calculations corresponds to the direct spin-phonon 
process, which is given by eq. (5). As previously explained, 
this process is screened out by PB. The calculated relax- 
ation rate follows the direct process until the temperature 
reaches about 0.8 K. Above this temperature, an exponen- 
tial growth of the rate is observed. This behavior would 
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Table 1: Levels involved in the Orbach mechanism for the two 
studied frequencies. The calculated strengths using DME as 
well as the fitted values to the data are summarized. 



not be expected for a pure two-level system. But within 
the giant spin approximation, the 5 = 4 Ni4 SMM has 
seven other levels lying above the two lowest ones, \a) and 
|6). Two of them, |a) and |/3), have been represented in 
fig. 2(a). Although they are weakly populated within the 
temperature range investigated, their presence can greatly 
enhance the relaxation rate. Instead of the direct relax- 
ation \b) — > |a), the system can be thermally excited from 
|6) to some high level, say \i), and then fall down to \a). 
This is the Orbach mechanism, yielding a characteristic 
Arrhcnius temperature dependence of the rate, i.e., the 
second term of cq. (6). Its activation energy Ei equals 
the separation between levels b and i, whereas its strength 
Foi depends on the coupling between \i) and \a) through 
transverse phonons. 

From the same statistical argument mentioned in the 
introduction, the Orbach mechanism can also be bottle- 
necked. But the emitted phonons have a much higher 
frequency Ei/h than those from the direct process, lead- 
ing to a much lower bottleneck parameter (i.e., there is a 
much higher density of phonon modes available for spin 
relaxation). Assuming the same inhomogeneous broaden- 
ing SiOo, Bi^g » 800 ioT Ei = 7 K, which is close to the 
energy separation between \b) and the fourth energy level 
for the two different applied fields. Assuming that Fph 
is independent of temperature, only processes with a rate 
larger than about 5 x lO'^ s~^ will be strongly bottlenecked 
(eq. (9)). This can be seen from the data, which behave 
differently at 10 and at 27.8 GHz in the high temperature 
region of fig. (4). Whereas the measured rate at 27.8 GHz 
is still two orders of magnitude lower than the calculated 
one, the data at 10 GHz tend to meet the DME calcula- 
tion. This is due to the fact that as the Orbach mechanism 
starts, the rate computed at 10 GHz using the universal 
DME is still low (F « 10^ - lO" s'^), and hence weakly 
bottlenecked. On the contrary, the relaxation at 27.8 GHz 
is still strongly bottlenecked in that temperature range be- 
cause the calculated rate is faster than 10* s~^. Table 1 
summarizes the three Orbach channels relevant for the en- 
ergy relaxation in the temperature range studied. It gives 
their activation energies, their calculated strengths using 
DME, and the values fitted to the data using eq. (6). The 
fitted values are smaller than the calculated ones because 
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PB is not completely suppressed. As explained, this effect 
is stronger at 27.8 GHz than at 10 GHz. 

It is interesting to estimate the phonon relaxation rate 
necessary to suppress the PB from eq. (9). If the phonon 
relaxation is very fast. B^g^T^s-[,,T ^ Tph, A_ = PspT, 
and there is no bottleneck. One would need to have Tph ^ 
10 ns at 10 GHz and Tph < 4 ns at 27.8 GHz to meet 
this requirement at low temperature. This would only be 
possible for a submicron-size crystal. 

In conclusion, the temperature dependence of the en- 
ergy relaxation rate between the two lowest-lying spin- 
states of Ni4 has been measured at two different frequen- 
cies. Using calculations of the rate based on the universal 
DME and taking into account PB that screens out the 
direct spin-phonon process at low temperature, the two 
sets of data are well accounted for in the whole tempera- 
ture range investigated. In particular, the Orbach mecha- 
nism dominates at high temperatures, and is weakly bot- 
tlenecked at 10 GHz. Finally, we point out that the fact 
that the universal DME calculations and the PB work so 
well to explain the data is surprising, as coherence is ex- 
pected to play an important role in collective spin-lattice 
relaxation in SMM crystals [15,39]. 
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